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Abstract 

Tensor models can be regarded as theories of dynamical fuzzy spaces, and pro- 
vide background independent theories of space. Their classical solutions corre- 
spond to classical background spaces, and the small fluctuations around them can 
be regarded as fluctuations of fields on them. In this paper, I numerically study a 
tensor model around its classical backgrounds of two-, three- and four-dimensional 
fuzzy flat tori and show that the properties of low-lying low-momentum modes are 
in clear agreement with general relativity. Numerical analysis also suggests that 
the lowest-order effective action is composed of curvature square terms, which is 
consistent with general relativity in view of the form of the considered action of 
the tensor model. 
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1 Introduction 



Tensor models (or generalized matrix models) jl]-[H] were originally considered as the general- 
ization of the matrix models, which describe the two-dimensional simplicial quantum gravity, 
to higher- dimensional cases. As in the matrix models, partition functions are given by a certain 
summation over simplicial manifolds, which are dual to the Feynman diagrams. Although the 
matrix models have successfully been analyzed extensively by powerful analytical methods, no 
analytical methods for tensor models are presently known. Moreover, partition functions of 
tensor models lack physical interpretations such as the topological sum over surfaces realized 
by the 1/iV-expansion in the matrix models. 

In ref. [9], a new interpretation of rank-three tensor models was proposed. This proposal 
is based on the facts that a fuzzy space is defined by an algebra of functions, and that an 
algebra can be characterized by a rank-three tensor that defines multiplication, f a fb = C a b c f c - 
In this new interpretation that a rank-three tensor may be regarded as the structure constant 
of an algebra defining a fuzzy space, rank-three tensor models can be regarded as dynamical 
theories of fuzzy spaces. This reinterpret at ion of the models is not just a formal difference, 
but provides a new practical way of extracting physics from tensor models. Namely, this 
new interpretation makes it meaningful to treat the models semiclassically, in which classical 
solutions of the models can be regarded as background fuzzy spaces, and small fluctuations 
around solutions as field fluctuations in fuzzy spaces. Another important difference from the 
original proposal p]-[8] is that three is enough as the rank of the tensor for describing a fuzzy 
space of any dimensions. This property simplifies the structures of tensor models, since various 
dimensional cases can be considered in one form. 

The rank-three tensor models have mainly been analyzed by numerical means by the 
present author, and several interesting results have been obtained [9J QUI dH [321 In 
particular, in the models that have a certain Gaussian type of classical solutions corresponding 
to fuzzy flat spaces of arbitrary dimensions, it has numerically been shown that the properties 
of low-lying low-momentum small fluctuation modes around the Gaussian solutions are in 
good agreement with the general relativity in the one-, two-, three- and four- dimensional 
cases [121 US] • In other words, general relativity was found to emerge in these models as an 
effective low-momentum description. However, in the three- and four- dimensional cases, owing 
to the heavy computational requirement, the agreement has been shown only for the numbers 
of the modes. To conclude the agreement, the main purpose of this paper is to carry out more 
detailed numerical analysis of the fluctuation modes in these dimensional cases by applying an 
approximate method similar to that in ref. [12J to the analysis of the tensor model introduced 
in ref. [13], which is much simpler than that analyzed in ref. |12j . 
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This paper is organized as follows. In the following section, the simpler model is reviewed. 
In £j3j by generalizing the discussions in ref. [13] to other dimensions, the metric fluctuation 
modes in general relativity are studied. In §H the results of numerical analysis for the two-, 
three- and four- dimensional fuzzy tori are shown and compared with general relativity. The 
final section is devoted to the conclusion and future problems. 



2 The model 



The tensor model that will be studied in this paper was introduced in my previous paper [T3] . 
I summarize the main points below. 

The dynamical variable is a totally symmetric rank-three real tensor: 

C abc : real, (1) 

C a bc = C bca = C cab = C bac = C acb = C cba . (2) 

There is also a nondynamical symmetric rank-two real tensor: 

g ab : real, (3) 
9 ab = 9 b \ (4) 

which is basically taken to be g ab = 5 ab (or 5 (a — b) for continuous indices). An action is 
assumed to have invariance under the orthogonal transformation M a b as 

C abc -> M a a 'M b b 'M c c 'C a , b , c , (5) 

This invariance can be achieved by contracting the upper and lower indices in pair, e.g., 
g aa> g bb ' g cc 'C a bcC a 'b>c>- In some places, I will use g ab (g ab = (g -1 )^) to raise (lower) the indices 
of tensors. 

The index of the tensors is the label of functions in a fuzzy space. It may take values in 
an infinite set or it may even become a continuous variable in describing fuzzy spaces with 
an infinite number of functions such as a fuzzy infinite flat plane, or may take values in a 
finite set for a compact fuzzy space. The former case will appear below in defining a tensor 
model which possesses strictly Gaussian solutions. The latter case will appear in where 
approximate Gaussian solutions corresponding to fuzzy flat tori are considered. The former 
case may be considered as an infinite limit of the latter case. 

At present, there are no principles for specifying an action. As an interesting case, I 
consider an action that has the Gaussian-type solutions 

C Xl x 2 x 3 = -Bexp [-/? ((xi - x 2 f + (x 2 - x 3 ) 2 + (x 3 - xi) 2 )] , (6) 
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g x ' X2 = S D (x 1 -x 2 ), 

where B and (3 are positive numerical constants, x« are D-dimensional continuous coordinates, 
Xi = (xj,x 2 , ■ ■ ■ ,xf), and (x) 2 = Y^=ii xf1 ) 2 - The summation in the contraction of pairwise 
indices is defined as f d D x. The algebra of functions f Xl f X2 — C XlX2 X3 f X3 defines a commutative 
nonassociative fuzzy flat space considered in ref. [H]. In the momentum basis, which is 
obtained by applying Fourier transforms to the indices, the solutions have the form 

C PlP2P3 =A5 D (pi + p 2 + p 3 ) exp [-a {p\ + p\ + p§) ] , (7) 
g™» = 6 D (p 1 +p 2 ) 1 

where A and a = are positive numerical constants. The summation in the contraction 
of pairwise indices is defined as J d D p. The main reason for considering such solutions is the 
(partial) computability due to the Gaussian forms. Another reason is that the fuzzy spaces 
are physically well-behaved in the senses that their fuzziness is well localized and that they 
have the Poincare symmetry. Moreover, as discussed in ref. [13J and as will be reviewed in 
the next section, one can consider a natural correspondence between metric fluctuations in 
general relativity and fluctuations in tensor models around the Gaussian solutions. 

There exist infinitely many actions that have the Gaussian solutions. An example was 
considered in ref. [12] , and the counting of low- lying low-momentum modes showed agreement 
with general relativity in the one-, two-, three- and four-dimensional cases. A much simpler 
action was considered in ref. [13] . This action contains a fractional power and cannot be 
considered to be well defined for arbitrary values of the tensors. However, at least for the 
small fluctuations around the Gaussian solutions, the action is well-defined and the properties 
of low-lying low-momentum modes showed clear agreement with general relativity in the two- 
dimensional case. In this paper, the analysis will be carried out for the much simpler action, 

S = W abc W abc , (8) 

where 

W abc = C abc - (iT-§)/(K-§)/(iH)/iW, (9) 

K a a ' = C ahc C a ' h \ (10) 

H a bc = C a deCbf d C c e f . (11) 



When A satisfies 



^=l^ f , (12) 



the Gaussian solutions ([7]) satisfy 

W abc = 0, (13) 
and become the classical solutions to the equation of motion from the action (jSJ). 
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3 General relativity on flat tori 



The strategy for comparing the tensor models with general relativity basically follows the 
procedure carried out for the two-dimensional case in my previous paper [IB] . Generalizing 
the Gaussian solutions ([6]) in an invariant manner, the fluctuations of metric tensor in general 
relativity and those of the tensor in tensor models are assumed to be related b}Q 

C X1X2XZ = Bg(xi)*g(x 2 )*g(x 3 )* exp [-(3 (d{x 1 ,x 2 ) 2 + d{x 2 ,x 3 ) 2 + d{x 3 ,x l ) 2 )] , (14) 
g XlX2 = S D (x 1 -x 2 ), 

where g(x) = det and d(x,y) denotes the distance between two points, x and y. 

Although it is possible to directly compare by the correspondence ( TT41) . which should give 
the same final results, the actual comparison will be carried out in terms of the tensor K ab 
defined in (fit)]) for technical simplicity, as in ref. [13]. In tensor models, a small fluctuation, 
SC a bc, around a classical solution, C° fcc , induces a fluctuation of K ab expressed as 

5K ab = 5C acd C\ cd + C° acd 5C b cd . (15) 



On the other hand, the correspondence (fT4l) induces 

6K pm = Sg^tj)! + p 2 ) (pi - p 2 T(pi - P2T ex P (~Y^( pi ~~ P2 ) 2 ) ' ( 16 ) 

where the momentum basis is used, and 5g^ v {p) is the Fourier transform of 5g liv {x)^ which is 
the metric fluctuation around a flat background. Equations ffTBI) and ffTB]) make it possible to 
compare the two distinct types of fluctuation, SC a bc in tensor models and Sg^x) in general 
relativity. 

In general relativity, not all of the components of g^x) are real geometric degrees of 
freedom on account of the gauge symmetry of local translations. As will be explained in the 
following section, the gauge symmetries in the tensor models appear as modes with vanishing 
spectra, and the modes to be compared are orthogonal to these gauge modes in the sense of 
the metric 

ds 2 c = 5C abc 5C abc , (17) 

which is invariant under the orthogonal transformation §5§ . As reported in my previous paper 
[13], by putting the correspondence f[T4")) into f|T7|) . one obtains the DeWitt supermetric [15], 



ds 2 DW = / d D x^gix) [(g^(x)5g^(x)f + Ag^(x)g^(x)5g tlv (x)5g pa (x)] . (18) 



Here g of tensor models should not be confused with the metric tensor g M „ (x) 
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Consequently, the metric fluctuation modes to be compared must be orthogonal to the local 
translational symmetry in the sense of the supermetric ([TBI) . 

In the following numerical analysis, D-dimensional fuzzy flat tori will be considered. There- 
fore, I will discuss below the orthogonal fluctuation modes in general relativity in the back- 
grounds of Z)-dimensional flat tori with metric = 5^. With such a background, the 
supermetric ffTS]) is explicitly given by 

D D 



ds^w — I d x 



[1=1 fi=X,v=2, 

fJ.<V 



(19) 



The infinitesimal gauge transformation in the flat background is given by 

8 v g^u{x) = <Vv(x) + dvV^x), (20) 

where v^(x) is the local translation vector, or 

S v g^(p) = iPnV v {p) + ipvV^(p) (21) 

in the momentum basis. 

As can be seen in ( l2TI) . at the vanishing momentum sector p^ = 0, the gauge symmetry 
is null, and all the components of the metric tensor are real geometric degrees of freedom. 
Diagonalizing the supermetric ( fl9l) . the modes are classified into three orthogonal classes: 

1 conformal mode: Sg^ = 5g uu ^ for all /i, v, others = 0, (22) 

D 

{D — 1) traceless diagonal modes: 5g^ = 0, off-diagonal components = 0, (23) 

fM=l 

D(D — l)/2 off-diagonal modes: diagonal components = 0. (24) 

At the nonvanishing momentum sector, one may take the momentum to be in the direction 
(pi, 0, • • • ,0) with evident generalization to the other directions. Then the explicit form of the 
infinitesimal gauge transformation (|2T|) is given by 

5 v g n = 2ip 1 v 1 , 

Svgifi = ipiVp for n ^ 1, (25) 
Kg^v = for /I, v ^ 1. 

In the same manner, the modes orthogonal to the gauge modes fl25l) in the sense of the 
supermetric ( fl9]) are classified into three orthogonal classes: 

1 diagonal mode: 5gn = —D + 1, 5g^ t = 5 for all /i ^ 1, (26) 



5 



D 

(D — 2) traceless diagonal modes: Sgn = 0, Sg^ = 0, off-diag. comps. = 0, (27) 
(D — 1)(D — 2)/2 off-diagonal modes: Sg^ ^ for fi ^ u, fi,v^l. (28) 

4 Numerical results 

I will first review the general strategy for obtaining the quadratic potential for small fluctua- 
tions in tensor models, which was carried out in refs. 



The components of C abc in the tensor models are subject to the symmetry constraint ([2]). 
The independent normalized fluctuation components are obtained from the metric f|T7|) as 



ds 2 c = * & > c )] SC abc 5C abc = SC (aM 5C^ c \ (29) 

(a,b,c) (a,b,c) 

where (a, b, c) denotes an order-independent set of a, b and c; C^^c) are the normalized inde- 
pendent components defined by C^( a ,6,c) — \/ m [( a ' ^ c )] ^a&cj and m[(a, 6, c)] is the multiplicity 
defined by 

{1 a = b = c, 
3 a = b ^ c, b = a, or c = a ^ b, (30) 
6 otherwise. 

Then the coefficient matrix of the quadratic potential for the normalized small fluctuation 
components around the Gaussian solutions ([7j) with (fT2"|) is given by 



]\/[(a,b,c),(d,eJ) _ 



dW ghi dW 9hi 



9C (aM) dC M) 



= Y \ dW 9hi dW9kl 

c=c y/m[(a,b, c)]m[(d,e,f)] dC abc dC def 



c=c 
(31) 

where I have used the form of the action (jSJ) and the fact that (1T31) is satisfied by the Gaussian 
solutions © with (fT2l) . 

The numerical analysis will be carried out around the backgrounds of D-dimensional fuzzy 
flat tori. These are the simplest fuzzy spaces with a finite number of functions on it and are 
suited for numerical and analytical analyses. Because of the translational invariance of the 
backgrounds, it is convenient to take momentum bases, where an index is assumed to take 
discrete values in a finite range as 

p=(p\ P 2 ,...,p D ), (p l = -L,-L + l, ■•■,£). (32) 

The Gaussian solutions ([7]) with evident replacement of the ^-function with a Kronecker 
delta will be the solutions to the action (jSj) in the limit L — > oo with 1/L 2 <a<l/L, since 
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the effects caused by the finiteness and discreteness of momenta will vanish in this limf||] 
In fact, in my previous paper [13], correct numerical solutions similar to the Gaussian solutions 
((7j) were explicitly found for the action ([8]) with such discrete and finite momentum indices in 
the two-dimensional case. Therefore, it can be expected that the coefficient matrix (13T1) with 
( !32|) will give a good approximation for D-dimensional fuzzy flat tori, when L is large enough 
and a is in the above range. 

The explicit form of the partial derivative in (l3"Tj) was given in [13] in the momentum basis 

as 



dW, 



abc 



m[(d,ej)] Y^"W(eK(/) 



(31)2 v(a)a{b)o(c) 
C=C cr,cr 



dCdef 

where the summations are over all the permutations of the indices, and 



1 a(a)a(b)a(c) ' \ OC> ) 



rrP4P5P(i —XP4.XP5XP6 _ OSP4+P5 XP6 ( I, U t, \~j r r 

P1P2P3 PI P2 P3 ° U P1+P2 U P3 \ n, Pl rv P2"'P3) U Pl,~P4,-pi+P4 U -ps,P2,-p2+P5 

-\~ 38p* P &p 2 +p 3 (kpi) 5 fpi,P2+P3 C Pl,~P4,~P5 C P2,P3,-P6 

+ 3(5^ S p2+P3 {k Pl ) s /p 1) p 2 +p 3 Cp 2j p 3) _p 2 _p3Cp 5 _|_p 6j _p 5) _p 6 , (34) 
c Pi,P2,p 3 =Aexp[-a(pl + p 2 2 + pj)}, (35) 

D 

k p =A 2 2 exp (Sap 2 ) , (36) 

{2 2 
fcpi -fc ?2 torfc Pl ^fc P2 , ^ (37) 
g(k Pl ) 9 for = /c P2 



with A given by (1121) . 

In general, the coefficient matrix (fSTl) contains a number of vanishing eigenvalues for a 
nontrivial classical background solution. These zero modes come from the broken generators 
of the orthogonal symmetry ([5]) in such a background, and have been analyzed completely 
for D-dimensional fuzzy flat tori in ref. [12]. Since these modes should physically correspond 
to local gauge symmetry nonlinearly realized in a symmetry-breaking background [16], they 
are not of interest in comparison with the orthogonal metric fluctuation modes discussed in 
the previous section. In the present approximate treatment, these modes generally have very 
small but nonvanishing eigenvalues and cannot be distinguished from the modes of interest 
particularly in the low-momentum regions, as was problematic in ref. [12]. To overcome the 
problem, in the present numerical analysis, the coefficient matrix (I3TT) will be projected to 
the subspace orthogonal to the approximate broken symmetry modes that are obtained by 



Tin this limit, the physical size of a torus becomes infinitely large in comparison with the scale of fuzziness 
it. 

■"■In the sense of the metric (flT)) . 
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Figure 1: Low- lying spectra on two-dimensional fuzzy torus for L = 10 and a = 1.5/L 2 . The 
horizontal axis is the size of momentum \p\ = y (p 1 ) 2 + (p 2 ) 2 of each momentum sector, and 
the vertical axis is the spectral value. The solid line is 1.4 x 10 _6 |p| 4 . 



infinitesimally applying the orthogonal transformations (J5]) to the Gaussian solutions ([7]). The 
eigenvalue /mode analysis will be carried out for the projected coefficient matrix. 

Because of the translational symmetry of the torus backgrounds, the coefficient matrix 
( l3Tj) is a direct sum of the submatrices of each momentum sector labeled by the momentum 
p = pi + P2 + P3 of the fluctuations 5C PlP2 p 3 . Therefore the eigenvalue/mode analysis will be 
carried out for each momentum sector of the fluctuations. 

The numerical analysis is performed on a Windows XP64 workstation with two Opteron 275 
(2.2 GHz, dual core each) processors and 8 GB memories. The C++ codes§ are compiled by the 
Intel C++ compiler 10.1 with OpenMP parallelization. ACML 4.0 is used for LAPACK/BLAS 
routines. Mathematica 6.0 is used for analyzing the outputs. 

Let me first show the results for the two-dimensional case. In my previous paper |13j . 

C+ 



^Interested readers can download the 
sasakura/codes/gausscodes.zip . 



codes from http://www2.yukawa.kyoto-u.ac.jp/' 
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Figure 2: Contour plots of 5K q - q+p for low-lying mode in p = (1, 0) sector in two-dimensional 
case. The left- and right-hand figures are (115j) from the numerical computations and (fl6l) 
substituted with ( 12"6l) . respectively. The axes are g = {q l ,q 2 )- 

an eigenvalue /mode analysis was carried out in a numerically correct way without using the 
approximation explained above. The merit of the approximation is that a larger L can be 
taken owing to the efficiency of computation. In Fig. HJ the case of L = 10 and a = 1.5/ L 2 is 
shown. In the figure, a trajectory of spectra can be clearly seen up to |p| < 6, and no overlap 
with some unwanted spectra at p — (2, 1) for L = 5 reported in my previous paper [13] is 
observed. This supports the expectation that, as L increases, the agreement with general 
relativity becomes clearer. The spectra are well located on the massless trajectory of \p\ 4 , 
which is in agreement with the result of [T3] . 

As reported in ref. [13], the properties of the low-lying modes are in clear agreement with 
the results from general relativity in the previous section. In the p = sector, there are three 
low-lying modes with the spectra 1.6 x 10~ 6 , 6.4 x 10~ 6 and 6.6 x 10~ 6 . By comparing the 
numerical results of (fT5|) and the analytical expression (TIB"]) , these modes were clearly identified 
with (12~4"1) . ( 1221) and (1231) . respectively. The low- lying mode at each nonvanishing momentum 
sector was clearly consistent with (T261) . as is shown for p = (1,0) in Fig. [2j 

In the three-dimensional case, I will show the numerical results for L = 3 and a = 1/L 2 . 
This number of L is the maximum value feasible within the machine memory size. In Fig. [31 
the low-lying spectra are plotted. \p\ A seems to be better than \p\ 2 as a fitting line for the 
spectra. 

In more detail, in the p = sector, there are three types of spectrum, 3.4 x 10~ 6 , 5.4 x 10~ 6 
and 7.7 x 10~ 6 with degeneracies of two, three, and one, respectively. Therefore, these modes 
can be consistently identified with the modes f[2"3"|) . (1241) and (1221 . respectively. In fact, the 
fitting of (fl6l) to (fl5l) of each mode from the numerical computations was clearly consistent with 
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Figure 3: Low-lying spectra on three-dimensional fuzzy torus for L = 3 and a = 1/L 2 . The 
solid and dashed lines are 8 x 10~ 5 |p| 4 and 5 x 10~ 4 |p| 2 , respectively. 




Figure 4: Contour plots of 5K q _ q+p for second low-lying mode in p = (1, 0, 0) sector in three- 
dimensional case. The left- and right-hand pairs of figures show (fl5l) from the numerical 
computations and (|T6l) substituted with (|26|) . respectively. The left- and right-hand figures in 
each pair show SK g _ q+p at the slices q = (q 1 , q 2 , 0) and q = (0, q 2 , q 3 ), respectively. 
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Figure 5: Low-lying spectra on four-dimensional fuzzy torus for L = 2 and a = 1/L 2 . The 
solid and dashed lines are 3 x lCT 4 |p| 4 and 1 x 10~ 3 |p| 2 , respectively. 



the identification. For example, fitting to the last mode resulted in Sg^ t = 2.45 ±0.02, Sg^ 
± 0.01, which is in clear agreement with 



In each p ^ sector, there are three low-lying modes. This number agrees with the result 
in the previous section. For example, the p = (1,0,0) sector has three low-lying spectra, 
1.4 x 10 -4 , 2.1 x 10~ 4 and 2.3 x 10~ 4 . By a similar fitting, they can be clearly identified with 
the modes (1281) . (|26|) and ff27|) . respectively. For example, an analysis of the second mode 
resulted in 5gn = —0.74 ± 0.01, 5(722,33 = 1-88 ± .01, 6g^ u = ± 0.01, which is in clear 
agreement with (12"oT) . In Fig. HJ f[T5"j) for the second mode and ([TBI) substituted with (12T)|) are 
shown. 

In the four-dimensional case, I will show the numerical results for L = 2 and a = 1/L 2 . 
This number of L is the maximum value feasible within the machine memory siz^l The 
low-lying spectra are plotted in Fig. [51 |p| 4 seems to be favored over \p\ 2 as a fitting line. 



^Even for L = 2, the precision of an array variable had to be made single to reduce the necessary memory 
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Figure 6: Contour plots of 5K q _ q+p for last low-lying mode in p = (1, 0, 0, 0) sector in four- 
dimensional case. The left- and right-hand pairs of figures show (fl5l) from the numerical 
computations and ( 1161) substituted with ( 1261) . respectively. The left- and right-hand figures in 
each pair show 5K q _ q+p at the slices q = (q 1 , q 2 , 0, 0) and q = (0, q 2 , q 3 , 0), respectively. 

In more detail, in the \p\ = sector, there are three types of spectrum, 1.1 x 10~ 6 , 
1.8 x 10~ 4 and 6.5 x 10~ 4 with degeneracies of six, one and three, respectively. Therefore, 
these can consistently be identified with the modes (|2"lj) . (I2"2"j) and (12"31 . respectively. In fact, 
the fitting of ( fl6l) to ( fl5l) from the numerical computations was clearly consistent with the 
identification. For example, an analysis of the second mode resulted in Sg^ = 0.683 ± 0.003 
and 5g pv = ± 0.003, which is in clear agreement with (1221) . 

In each \p\ ^ sector, there are six low-lying spectra including their degeneracies, which 
agrees with the result in the previous section. For example, in the p = (1, 0, 0, 0) sector, 
there are three types of spectrum, 2.4 x 10~ 4 , 3.5 x 10~ 4 and 5.1 x 10~ 4 with degeneracies 
of two, three and one, respectively. Therefore these can be identified with the modes (I27p . 
(I28j) and (1261) . respectively. In fact, similar fitting procedures for each mode were consistent 
with this identification. For example, an analysis of the last mode resulted in 5gn = —0.41 ± 
0.01, £#22,33,44 = 0.80 ± 0.01 and = ± 0.01, which can be well identified with ([26]). In 
Fig. El (|T5l) of the last mode and ( Till substituted with ( 1261) are shown. 

In summary, the numbers and properties of the low-lying modes of tensor models on fuzzy 
flat tori are numerically shown to be fully consistent with general relativity, and the momentum 
dependence of the spectra is found to be \p\ 4 . The latter fact suggests that the lowest-order 
effective action is composed of curvature square terms. 

5 Conclusion and future problems 

In this paper, I have considered a tensor model that possesses Gaussian classical solutions 
corresponding to fuzzy flat spaces of arbitrary dimensions. The numerical analysis of the small 
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fluctuations for two-, three- and four- dimensional fuzzy flat tori have shown clear agreement 
with general relativity on these tori. Therefore, it can be concluded that general relativity 
appears as an emergent phenomenon in this model. 

The tensor models have a rank-three tensor as their only dynamical variable, and they 
are formulated background independently. Various background spaces can be generated as 
classical solutions 0, QUI E], and the above conclusion shows that small fluctuations around 
them can be identified with general relativity at least in a certain class of tensor models. Thus, 
the tensor models can provide an interesting unifying description of space and geometry. The 
specific choices of actions studied thus far must be generalized in future studies. 

In the studies thus far, the actions were assumed to have certain semipositive definite 
forms and the classical solutions were assumed to be their absolute minima. However, these 
assumptions force the effective actions to start with curvature square terms as is consistent 
with the numerical analysis presented in this paper, since the Euclidean Einstein-Hilbert 
action is topological in two dimensions and generally takes both positive and negative values 
in dimensions higher than two [17J. Therefore, indefinite actions or stationary solutions will 
be more interesting for the realization of the Einstein-Hilbert action, and will also provide 
further consistency checks in view of the qualitative difference between the two- and higher- 
dimensional cases. 

The next interesting question would be which fields can emerge in tensor models. Prob- 
ably, the gauge fields can emerge from an emergent metric field by a type of fuzzy analogue 
of the Kaluza-Klein mechanism. Then how about fermions? In fact, tensor models have an 
interesting potential. In the analyses performed thus far, the rank-three tensor is restricted to 
be totally symmetric, which requires functions in fuzzy spaces to be commutative. By simply 
relaxing this restriction to another milder condition such as the generalized Hermiticity con- 
dition stated in refs. [HO [3], the functions can also become anticommutative, i.e., fermionic. 
These ideas suggest some interesting directions of future study in this area. 
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